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We re-examine the constraints on the cosmic string tension from Cosmic Microwave Background 
(CMB) and matter power spectra, and also from limits on a stochastic background of gravitational 
waves provided by pulsar timing. We discuss the diflerent approaches to modeling string evolution 
and radiation. In particular, we show that the unconnected segment model can describe CMB 
spectra expected from thin string (Nambu) and field theory (Abelian-Higgs) simulations using the 
computed values for the correlation length, rms string velocity and small-scale structure relevant 
("^ ■ to each variety of simulation. Applying the computed spectra in a fit to CMB and SDSS data we 

^S] ' find that Gfi/c^ < 2.6 x 10^^ (2cr) if the Nambu simulations are correct and G^i/c^ < 6.4 x 10~^ 

in the Abelian-Higgs case. The degeneracy between G^/c? and the power spectrum slope ng is 
^'. substantially reduced from previous work. Inclusion of constraints on the baryon density from 

Big Bang Nucleosynthesis (BBN) imply that ns < 1 at around the Aa level for both the Nambu 
and Abelian-Higgs cases. As a by-product of our results, we find there is "moderate-to-strong" 
Bayesian evidence that the Harrison-Zel'dovich spectrum is excluded (odds ratio of ~ 100 : 1) by 
the combination of CMB, SDSS and BBN when compared to the standard 6 parameter fit. Using 
the contribution to the gravitational wave background from radiation era loops as a conservative 
r^ , lower bound on the signal for specific values of G^,/(? and loop production size, a, we find that 

^ . G^i./c^ < 7 X lO"'' for (xc^/(VG^i) < 1 and G^i/c^ < 5 x 10-"/a for ac^/{VGii) > 1. 

r^ ' PACS numbers: 98.80.Cq, 98.80.Jk 

o 

^ . I. INTRODUCTION 

The idea that cosmic strings may have been a source of inhomogeneities in the Universe has been around for 
V^J \ some time (see e.g. ref. [l| for an overview). Even though it now seems that inflation is the primary source of the 
inhomogeneities which lead to anisotropies in the Cosmic Microwave Background (CMB), the idea has persisted. The 
p^' ' main reason for this is that the amplitude of fluctuations produced by strings formed close to a Grand Unified Theory 
t:;^ , (GUT) phase transition {^ lO^^GeV) is of the order of observed density inhomogeneities on large scales, making it 
^) • an attractive probe of physics at energies way above that possible in terrestrial accelerators. 

It may seem very non-minimal to have two mechanisms creating inhomogeneities and contrary to the common-sense 
rS ' approach advocated by Occam's Razor. However, there are at least two classes of inflation models, linked to popular 
f^ I ideas in fundamental physics, where strings can be formed and the two sets of perturbations have similar amplitudes. 
One is based on Supersymmetric Hybrid Inflation [2-Q and the other is Brane Inflation [3]. In both cases cosmic 
strings can be formed during the phase transition at the end of inflation [8l-[l(j| . It was shown in ref. Ill that this can 



happen in a large number of symmetry breaking schemes and it is possible that semi-local strings could form J12l . Il^ 
in situations where strings arc not forced to exist on topological grounds. 

One cannot rule out strings using observations per se, rather one can only constrain the tension (or the mass 
per unit length) of the strings, and hence the energy scale at which the strings arc formed. The tension is usually 
quantified in terms of the dimensionless ratio Gfx/c^. This can be done in many different ways. In this paper we will 
focus on two: the anisotropies, created by the Kaiser- Stebbins effect [ij, [l^ and related pre-recombination processes 
in the CMB; and the stochastic background of gravitational waves created by the decay of string loops piij . 

There are serious difficulties in setting such limits due to the dynamic range involved in modeling the evolution of 
a cosmic string network over cosmological scales. This has led to a number of different constraints being published 
in the literature which are often contradictory - these will be reviewed in the subsequent sections. The objective 
of the present paper is two-fold. For both types of observations we will review the important physics and attempt 
to delineate the key uncertainties. We then present what we believe are the most reliable constraints based on the 
present knowledge. In the case of the CMB observations we will present two different constraints, based on different 
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FIG. 1: (Left) Comparison of cosmic string power spectra computed with the Unconnected Segment Model (USM) using the 
Nambu - model A - (solid line) and AH mimic - model B - (dotted line) parameters. Spectra have been normalized to the 
WMAP value of Cio, giving G/i/10~^c^ = 1.18 for the Nambu model and 1.91 for the AH mimic. (Right) Comparison of scalar, 
vector and tensor modes for the USM AH mimic (solid) parameters and the actual AH spectra from simulations (dotted) [25l |. 
The former has G/i/lO^'^c^ — 1.91, and the latter 2.04. From top-to-bottom at low £ the ordering of the curves is the total 
power, then the anisotropy from vectors, scalars and tensors, respectively. 

assumptions about the evolution of the string network. We note that the constraints presented here will apply only 
directly to strings where the intercommutation probability, p, is unity; that is, they can only be directly applied to 
normal field theoretic strings and not necessarily to cosmic supcrstrings [l%. 



II. CMB SPECTRA FROM COSMIC STRINGS 



The key to predicting the observed CMB anisotropics from strings is a detailed understanding of the evolution of a 
cosmic string network. The main idea is that, after formation, the network evolves toward a self-similar scaling regime, 
whereby the average properties of the network remain approximately constant as a function of time. In the standard 
scenario this is achieved by the production of loops and their subsequent decay into radiation, usually assumed to be 
gravitational. 

The scaling assumption is not at issue. It has been observed in just about every simulation that has ever been 
done. However, the precise characteristics of the scaling regime are of critical importance for predicting the CMB 
anisotropy spectrum. The problem is that each simulation technique needs to make assumptions in order to achieve 
sufficient dynamical range. The two most popular approaches are: (A) to solve the Nambu equations of motion for 
a string in an expanding universe and include the effects of reconnection by hand; and (B) to solve the equations 
of motion for the Abclian-Higgs (AH) model. In method A one ignores the effects of radiation backreaction, which 
is likely to be important in setting the size of loops which the string network creates, making the observed scaling 
dependent on the removal of sufficiently small loops from the simulation. Radiation into propagating modes of the 
field is included in method B, but one is necessarily forced to consider simulation box sizes which are only '^ 300 
times greater than the core string width, thereby reducing the dynamical range. This could also lead to the radiation 
being stronger in the simulation than in realistic string networks since the typical radius of curvature is large relative 
to the core width. Moreover one has to accept contortions involving fixing, or reducing the rate of growth of, the 
core width in comoving units in order to perform simulations in an expanding universe. Dynamical range is less of a 
problem in method A, since the core width is zero and the strings arc represented as ID objects, as opposed to the 
3D box required for field theory simulations, thereby allowing much larger grid sizes. Having said that no simulation 
can achieve the dynamical range required to model the real Universe! 

Simulations using method A have been used to derive the key statistical parameters of the string network which 
appear to be different in the radiation and matter eras |18l - (21| . It is found that the correlation length L = ^dii{t), where 
duit) is the particle horizon distance, differs substantially between the radiation and matter eras, with ^rad = 0.13 and 
^mat = 0.21 (see refs . [22,[23| and references therein), corresponding to a significantly lower string density (pstr — t^/L^ 
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FIG. 2: A montage of plots which show the dependence of the string spectrum on parameters of the USM. In each case the solid 
and dotted lines are the Nambu and AH mimic models respectively from Fig. [T] In all but the top-left plot we have varied the 
one parameter keeping all the others the same as in the AH mimic model: (top-left) the dashed line is the USM spectrum with 
^ — 0.13, V = 0.65c and /? — 1.9 but no matter-radiation transition; (top-right) the short dashed line has ^ = 0.2 and the long 
dashed line has ^ — 0.5; (bottom-left) the short dashed line has v = 0.0c and the long dashed line has v — 0.8c; (bottom-right) 
the short dashed line has 13—1.5 and the long dashed line /3 = 1.9. 



where /i is the string tension - whose value is the main subject of the paper) in the matter era. There is substantial 
small-scale structure on the long string network. This can be quantified in terms of a renormalized string tension |2J] , 
fJ-eS = /?M, and it is found that /3rad = 1-9 and /3mat = 1-5. On the other hand, the RMS velocity of the string segments 
only varies very slightly between the different eras, with Vrad = (''^rad)^^^ = 0.65c and Wmat = ("^^mat)^^^ = 0.60c. It has 
been shown that a velocity dependent one-scale model can be used to describe the evolution between the radiation 
and matter eras, and the extrapolation into the cosmological constant dominated era. This typically involves a long 
transient with the string density only achieving its matter era value around z ~ 100. 

Contrasting results arc found using method B [25|. Typically very small loops are created from the fragmentation 
of horizon scale loops and most of the energy loss, required in order to maintain scaling, is achieved by radiation of 
the massive modes associated with the AH model, f « 0.3 and v w 0.5c [26| are computed from the simulations, with 
little difference observed between simulations performed in the radiation and matter eras. The nature of the observed 
small-scale structure on the strings is different but is compatible with /3 w 1. 

Wc note that Ringeval et al [231 ^^'^ ^^^^^ dad ~ 0.16 and ^mat ~ 0.19 using Nambu simulations which corresponds 
to a change in density between the radiation and matter densities of « 1.4 as opposed to ~ 2.6 which would be the 
case if ^,ad ~ 0.13 and ^,nat ~ 0.21. This will lead to results more similar to those which come from the AH model. 
In addition an attempt to compute the spectrum directly from Nambu simulations has recently being made j28| . 

In both cases one can in principal compute the unequal-time correlator (UETC), C^^ypa = {dtj,u{T)Ofj^i,{T')), of the 
string network stress-energy tensor, 0^1^, from the simulation. Using the stiff source approximation, one can then 
compute the CMB anisotropics by either decomposing the UETC in eigenmodes [231, or averaging realizations with 



the correct statistical properties. An alternative approach is to model the UETC using what has become known at 
the unconnected segment model (USM) . This models the string network as an ensemble of string segments of length 
^du(t) with an RMS velocity v. This idea [30| was first used in the context of the CMB in refs. [311,1321 and was later 
adapted to include the matter-radiation transition [22| and the effects of small-scale structure on the network via the 
'wiggliness' parameter, /? [23|. The weakness of this approach is that there is some freedom due to the parameters of 
the model, ^, v and /3, which could be functions of time, but these can be and have been computed from simulations. 

The present work uses the USM and the code CMB ACT [23 which was also used in refs. [33l - [37| . The spectrum 
computed using the USM and the parameters from the Nambu simulations (model A) are presented in the left hand 
panel of Fig. [1] along with the spectrum computed from the AH model simulations [23 ■ It is clear that the two spectra 
are very different, with model A having a fairly sharp peak at around £ ~ 500, whereas the AH spectrum has a lower 
and broader peak, and falls off very sharply at high i. It is not necessary for the string spectrum to be known to the 
level of precision of the inflationary component, since it is constrained to be < 10% of the total anisotropy. However, 
the discrepancy illustrated here is larger than is acceptable and it leads to important quantitative differences when 
each of the spectra are applied to the data. 

The spectrum computed using the USM is not unique, being a function of the parameters, and in particular the 



first CMB spectrum computed using it j31| looks very similar to that computed from the Abelian-Higgs simulations. 
This first spectrum just used a single, constant value for ^ and u, and ignored the effects of the matter-radiation 
transition and small-scale structure, that is, it had (3=1. It seems plausible that it might be possible to compute an 
AH 'mimic' model using the USM (model B). In order to investigate this we ran a suite of simulations of the USM 
with 3 parameters, £,, v and (3, all assumed to be constant with time, and then minimized the residuals to the true 
AH spectrum taking into account the decomposition into the scalar, vector and tensor (SVT) modes, all of which are 
present in cosmic string spectra. We found that ^ = 0.35, v = 0.4c and (3 = 1.05 gave the best fit. These are quite 
close to the values computed from the AH model simulations discussed earlier. The spectrum created using model B 
is compared to that from the AH model simulations in the right hand panel of Fig. [T] We see that there is remarkably 
good agreement between the spectra and that the values of G^/c^ required to achieve the observed level of fluctuations 
are very similar. This suggests that the discrepancy in the spectra between the USM used to describe the Nambu 
simulations and that from the AH model simulations is mainly due to the computed correlation properties, and the 
differences in the simulation techniques, not the subsequent methods used to infer the CMB spectrum. We have also 
computed the CDM power spectrum for both string models, but these are subdominant compared to the inflationary 
component when normalized to the observed large angle CMB amplitude ,2%, i3l| . Since the string contribution is 
limited to < 10% of the total anisotropy, the CDM power spectrum from strings is completely subdominant. 

It is worth investigating what are the primary reasons for the differences between the Nambu and AH spectra. In 
Fig. [2] we have first removed the matter-radiation transition, making the spectrum just a function of single values of ^, 
V and (3. We have then varied each of the three parameters relative to the AH mimic model. It is clear that the most 
significant reason for the discrepancy is the effect of the matter-radiation transition, which is present in the Nambu 
simulations but not observed in the AH simulations. The other important effect comes from the correlation length 
^, whereas the parameters v and (3 appear to have little effect on the overall spectrum. The small-scale structure 
parameter, /3, docs play an important role in regulating the relative contribution of the scalar, vector and tensor 
components; values of /3 > 1 have an enhanced scalar component relative to the vector component. The tensor 
component is low in all cases. 

In the next section, where we describe constraints on Gfi/c^, we will present results for the two spectra (model A 
and B) - the qualitative conclusions are similar, but there are important quantitative differences, particularly in the 
upper limit on Gixjf? . It is instructive to try understand why there is such a significant discrepancy in the properties 
derived from the 2 different simulation methods which leads to these differences. 

The most worrying aspect of the AH simulations is their size relative to the core width. It could be that the absence 
of small-scale structure is related to the artificially strong massive radiation seen in the simulations. It is plausible 
that this could smooth out the small-scale structure on the long strings and remove the need for loop formation. The 
fact that the string width is kept constant in comoving coordinates, which makes it artificially large relative to the 
horizon, may be responsible for there being little difference between the scaling behaviours in the radiation and matter 
eras. The achilles heel of the Nambu simulations is the absence of radiation which could over-produce small-scale 
structure. Clearly more work needs to be done to resolve these discrepancies, but it is reassuring that both models 
can be represented with the necessary level of accuracy using the USM. An alternative to presenting results for the 
two separate models would be to allow ^, v and f3 to vary within the USM and marginalized over them as nuisance 
parameters. 





Model 


Parameter 


NS 


S (NAMBU) 


S (AH) 


HZNS 


HZ S (NAMBU) 


HZ S (AH) 


fib/i" 


0.0231 ± 0.0005 


0.0239 ± 0.0008 


0.0240 ± 0.0008 


0.0243 ± 0.0004 


0.0253 ± 0.0007 


0.0253 ± 0.0005 


Q,h^ 


0.107 ±0.003 


0.107 ±0.004 


0.107 ±0.004 


0.105 ±0.003 


0.105 ±0.003 


0.105 ±0.003 


Ba 


1.042 ± 0.002 


1.043 ± 0.002 


1.044 ± 0.003 


1.046 ± 0.002 


1.047 ± 0.002 


1.047 ± 0.002 


TR 


0.087 ±0.017 


0.088 ±0.017 


0.089 ±0.017 


0.113 ±0.018 


0.111 ±0.018 


0.101 ±0.017 


log(10i»P,j) 


3.04 ± 0.04 


3.01 ±0.04 


3.02 ± 0.04 


3.11 ±0.04 


3.07 ± 0.04 


3.04 ± 0.04 


Us 


0.957 ±0.013 


0.958 ±0.013 


0.970 ±0.016 


1 


1 


1 


Gm/IO-^c^ 
(/lo) 


- 


<2.6 

(< 0.044) 


<6.4 
(< 0.093) 


- 


1.9 ±0.5 
(0.027 ±0.013) 


5.3 ±1.0 
(0.075 ± 0.025) 


/i 


0.74 ± 0.02 


0.75 ±0.02 


0.75 ±0.02 


0.77 ±0.01 


0.79 ± 0.02 


0.78 ± 0.02 


-log£ 


1440.8 


1440.2 


1440.3 


1446.7 


1445.2 


1442.2 


+BBN 


1442.0 


1442.0 


1442.1 


1450.4 


1450.5 


1448.4 



TABLE I: CMB ± SDSS constraints on string models using gatr as the parameter. Models denoted by S use strings and NS 
indicates that we have not used strings in the MCMC, while HZ denotes a Harrison-Zel'dovich spectrum (i.e. ris = 1). NAMBU 
refers to the USM spectrum using the parameters derived from the Nambu simulations and AH to that of the AH mimic model. 
Upper limits on G^/c? and /lo are 2cr. The x^ values of the models can be computed from x^ ~ — 21og£. 





CMB 


±SDSS 


±BBN 


±SDSS/BBN 


NAMBU 


2.8 (0.054) 


2.6 (0.044) 


2.2 (0.032) 


2.2 (0.030) 


AH 


6.8 (0.11) 


6.4 (0.093) 


5.1 (0.057) 


5.0 (0.055) 



TABLE II: 2(T upper limits on the string tension G/x/lO "^ (? and /lo (bracketed) using gstr as the parameter. 



III. CONSTRAINTS ON G^jl/c^ FROM CMB ANISOTROPIES 



A. Methodology 



Our methodology for obtaining CMB constraints on G^jl/c^ will follow our previous work j36| which focussed on 
the implications of the WMAP 3 year data release. Since the temperature power spectrum produced by strings is 
limited to around ~ 10% of the total observed power, it is sufficient to keep the shape of the string spectrum fixed in 
a Markov-Chain- Monte-Carlo (MCMC) analysis, and vary only the normalization G^/(? . In this way, one avoids the 
time-consuming computation of the string spectrum for each set of cosmological parameters. To generate the fixed 
string spectrum for the Nambu and AH mimic (hereafter we refer to this simply as 'AH') models outlined above, we 
used best-fit parameters from the WMAP 5 year data release 38| . 

To generate the MCMC chains we used COSMOMC [H, which makes use of the CAMB [H and CMBFAST [U 
codes. We vary a total of 7 cosmological parameters: - the baryon density, Jlb/i^; cold dark matter density, Vtch?\ 
the acoustic scale. Ox] optical depth, r; amplitude and spectral index of primordial fluctuations, Pr and ns; together 
with the string normalization G^/c^. We will include Gjjl/c^ in two ways. In one set of runs we adopt Qstr = 
(G/i/1.1 X 10~^c^)^ as the parameter. This was used in ref. [i^l, and is similar to the method of ref. [4j], which used 
/lo, defined as the fraction of total power due to strings at £ = 10. An alternative parameterization is log]^o(^/^/"^^): 
which was used in ref. [36| . As we will see there are subtle differences in the resulting constraints because the choices 
impose flat priors on different functions of G^/c^. In both approaches we will compute /lo for each sample in the 
Markov-Chain in order to compare with other analyses of st ring constraints |43|. Finally, we marginalize over the SZ 
amplitude ^sz, assuming the Komatsu and Seljak template |44{ . Also, in some instances we will fix ng = 1, that is a 
Harrison-Zel'dovich (HZ) spectrum, to compare the goodness-of-fit to a model where ng is allowed to vary. 

We use CMB data from the WMAP 5 year release gl 13 and the ACBAR [13, BOOMERANG 0, CBI [13 and 
QUAD [4^ experiments, which observe to higher L The latter high resolution data arc important to help break the 
degeneracy which exists between ns and GpLJc? for WMAP alone. Following ref. [43, we also include the BIMA [50| 
data point at £cff ~ 5000 to constrain the SZ amplitude, without assuming a prior upper limit on ^sz as done in the 
analysis of the WMAP 5 year data [4a]. However, imposing this upper limit on ^sz docs not significantly effect our 
conclusions. We also use matter power spectrum data from the SDSS Luminous Red Galaxies sample [5l|, and Big 
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FIG. 3: 2D likelihoods of scalar spectral index ns versus cosmic string amplitude parameter gatr for the Nambu (red, likelihood 
surface at front) and AH (blue, surface behind) string spectra. In the left panel we show constraints from CMB+SDSS data, 
and on the right constraints from CMB+SDSS+BBN. 

Bang Nucleosynthesis (BBN) constraints on the baryon fraction from measurements of deuterium at high redshift 53 . 
The latter gives Vlhh^ ~ (2.13 ± 0.10) x 10~^, which is slightly lower than that suggested from CMB data alone 38| . 



B. Using QstT as a parameter 



The detailed results of our analyses using CMB+SDSS data and (/str as the parameter describing the amplitude of 
the string spectrum are given in table HI We include values for — log£ both with and without the BBN prior. We 
will use the results from CMB+SDSS data as our headline results in this paper, but will also discuss the significant 
impact of including BBN data. 

For models where ns is allowed to vary, the improvement in the goodness-of-fit when including strings is not 
significant over the non-string model, in both the Nambu and AH cases (Ax^ = 1.2 and 1.0 respectively for one extra 
parameter in the case where we do not include the BBN prior). This suggests that there is no need to include the 
strings in order to explain the data, but constraints can be imposed on the cosmic string tension. 

The constraints on Gji/c^ and /lo are summarized for different combinations of data and the two string models in 
table nil It is clear that constraint is improved by the inclusion of the SDSS and BBN data. The 2a upper limits 
which wc will quote from CMB+SDSS data arc G^/c^ < 2.6 x 10~^ for the Nambu model and G^/c^ < 6.4 x 10^^ 
for the AH model. The upper limit on G^l/c^ (and the respective string fraction /lo) is around a factor of two and 
a half smaller in the Nambu case and there are two reasons for this. The most significant is that the Nambu model 
spectrum has a lower value of Gp^j (? for a given power level on large-scales. In addition the Nambu spectrum has 
more small scale power relative to the AH case. This means that, if the amplitude of the spectrum is too high, the 
strings can modify the precise peak structure in the high WMAP signal-to-noise regime of £ « 200 — 400, leading to 
a tighter constraint. 

The degeneracy between G^jl/c^ and ng is illustrated in the left panel of Fig. |31 In previous work, wc found 
that ng = 1 was completely compatible with CMB data in the presence of a subdominant string contribution [36| ; 
this point was also made in ref. [43. With improved CMB data and the inclusion of the SDSS, we find this is no 
longer the case, especially for the Nambu model. There is still a degeneracy between Gp and ns in the AH model, 
although it is considerably reduced, but the two quantities appear to be completely decoupled in the Nambu case. 
The marginalized value of ng is around ia from ng = 1, with Ax^ = 10.0 between the HZ spectrum and variable 
ng models (for comparison, with no strings, Ax^ = 11.9). For the AH case, the marginalized value of ng is around 
2cr from ns = 1, with Ax^ = 3.8. This point has led us to suggest that the F- and D-term supersymmetric hybrid 
inflation models are very tightly constrained and that minimal models are probably ruled out by the data even in the 
presence of strings [Sj] . 





Model 


Parameter 


NS 


S (NAMBU) 


S (AH) 


HZNS 


HZ S (NAMBU) 


HZ S (AH) 


fib/i" 


0.0231 ± 0.0005 


0.0234 ± 0.0007 


0.0234 ± 0.0007 


0.0243 ± 0.0004 


0.0248 ± 0.0007 


0.0251 ±0.0005 


Q,h^ 


0.107 ±0.003 


0.107 ±0.004 


0.107 ±0.003 


0.105 ±0.003 


0.105 ±0.003 


0.105 ±0.003 


Ba 


1.042 ± 0.002 


1.043 ± 0.002 


1.043 ± 0.002 


1.046 ± 0.002 


1.046 ± 0.002 


1.046 ± 0.002 


TR 


0.087 ±0.017 


0.087 ±0.016 


0.088 ±0.017 


0.113 ±0.018 


0.112 ±0.019 


0.103 ±0.018 


log(10i»P,j) 


3.04 ± 0.04 


3.02 ±0.04 


3.03 ± 0.04 


3.11 ±0.04 


3.09 ± 0.04 


3.05 ± 0.04 


Us 


0.957 ±0.013 


0.957 ±0.013 


0.960 ±0.014 


1 


1 


1 


Gm/IO-^c^ 
(/lo) 


- 


<2.3 
(< 0.035) 


<5.3 

(< 0.071) 


- 


<2.5 
(< 0.044) 


(4.7 ±1.4) 
(0.062 ± 0.028) 


/i 


0.74 ± 0.02 


0.74 ±0.02 


0.74 ± 0.02 


0.77 ±0.01 


0.79 ± 0.02 


0.78 ± 0.02 


-log£ 


1440.8 


1440.2 


1440.4 


1446.7 


1445.2 


1442.3 


+BBN 


1442.0 


1442.1 


1442.0 


1450.4 


1450.6 


1448.5 



TABLE HI: Equivalent of tableUbut using logj^Q(G/x/c^) as the parameter defining tlie string amplitude. 





CMB 


±SDSS 


±BBN 


±SDSS/BBN 


NAMBU 


2.5 (0.042) 


2.3 (0.035) 


1.8 (0.021) 


1.8 (0.020) 


AH 


5.7 (0.084) 


5.3 (0.071) 


3.9 (0.038) 


3.8 (0.035) 



TABLE IV: Equivalent of table HIl but using logjQ (G/i/c^) as the parameter defining the string amplitude. 



The values of /lo we obtain from our analysis constrain the string component to be less than 4.4% of the large scale 
amplitude in the case of the Nambu model and 9.3% for the AH model, using the combination of CMB + SDSS data. 
In the AH case the values we find are smaller than those quoted in ref. [4^]. For CMB only, they found a nominal 
2(7 preference for strings with /lo = 0.11 ± 0.05, whereas we only find a 2cr upper limit of 0.11, corresponding to 
G\il(? < 6.8 X 10~^. We presume that this is as a result of the improved CMB data. 

It is noticeable that imposing the BBN prior reduces the upper limit on Gfi/c^ and /iq. This is because higher 
values of G^l/c^ are degenerate with larger values of i7b/i^ and ns, as illustrated in ref. [S^- This effect is more 
pronounced when using the AH model where there is still some degeneracy without the inclusion of the BBN prior. 
Applying the prior entirely removes the degeneracy between G^/c^ and ng, as illustrated in the right panel of Fig. [31 
In this case, for both the Nambu and AH mimic models we find ng < 1 at around 4ct significance (Ax^ = 16.8 for the 
Nambu model and 12.8 for the AH model). 



C. Using logio(G/i/c^) as a parameter 



For completeness we have included parameters and limits derived from using \ogif^{G fi / c^) as the parameter to 
describe the amplitude of the string power spectrum. The results are summarized in tables IIIII and IIVI along with 
Fig. m The results are broadly consistent with those obtained using gstr as the parameter, but there are some 
quantitative differences which come from the different prior space. This means that there is a different measure used 
in the marginalization process. We note that log]^g(G'/.t/c^) was used in our earlier work [36|, whereas something closer 
to gstr was used in other previous work and this accounts for some of the discrepancies with the results reported in 
ref. [43|. There is, of course, no correct choice of parameter and therefore the differences represent the level of "natural 
uncertainty" in the parameter estimation process. 

In this case we find that the upper limits are modified to Gy^jt? < 2.3 x 10'^ for the Nambu model and < 5.3 x 10"^ 
for the AH model using CMB±SDSS, and limits from other data combinations are also strengthened. There are also 
some small differences in the marginalized values of cosmological parameters which is a result of the different weighting 
in the marginalization process. These differences are slightly larger for the AH model where there is more significant 
degeneracy between the cosmological parameters and Gy/c^. As expected, the values of — log£ are more or less 
unchanged between using ^str and logjo(G/.t/c^). 
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FIG. 4: 2D likelihoods of ns versus logj^jJG/^/c^]. Labeling is the same as in Fig. [S] 



D. Bayesian evidence against the Harrison-Zel'dovich spectrum 



As a by-product of our analysis, we find there is strong evidence that ng 7^ 1 within the 6 parameter family of 
ACDM models. With no strings, including the BBN prior increases the Ax^ between HZ and variable rtg from 11.9 to 
16.8 (as noted above, a similar Ax^ for BBN is also true with strings). One can further quantify this significance by 
computing the Bayesian evidence using the Savage-Dickey method (see, for example, ref. |5J]). To do this, we assume 
Gaussian prior information on ng with mean 1 and standard deviation 0.2. The resulting Bayes factor between HZ and 
variable ns models is logi? = —4.5, corresponding to an odds ratio of ~ 100 : 1, which provides "moderate-to-strong" 
evidence that ?7g ^ 1. 



IV. CONSTRAINTS ON Gfi/c^ FROM PULSAR TIMING 



It is expected that the cosmic string network will emit gravitational waves, the dominant contribution coming 
from loop decay [T^l, leading to a stochastic background which can be constrained by pulsar timing and also recent 
measurements made by LIGO |56l - [63 |. At this stage those from pulsar timing are the strongest, and therefore we shall 
concentrate on them. Just as with the limits from the CMB anisotropy there are a number of constraints available in 
the literature, many of which are contradictory, mainly due to the different assumptions which are made. 

There are two important issues here: the size of the loops produced by the string network and the spectrum of the 
radiation they produce. These issues require detailed modeling of the behaviour of the string network on small scales 
and an understanding of the impact of radiation backreaction. 

The size of loops produced is usually assumed to be a constant fraction of the horizon size at formation, a = 
l{t{)/d-H{t), and will be related to the gravitational backreaction scale by a '^ F Gfj,, where F '^ 60 is defined below. The 
Nambu simulations do not include the effects of radiation backreaction, and we have argued that the AH simulations 
include massive radiation at an unnaturally strong level. Therefore, there is no reliable determination of the loop 
production size a from either type of simulation. 

The spectrum of radiation emitted by given loop solutions of the Nambu equations of motion can be computed using 
the quadrupolc formula under the assumption that backreaction does not effect the shape of the string trajectory. 
In particular one finds that the total power emitted by a loop is independent of the size the loop, I, and if given by 
P = T G^^c. This can be written in terms of harmonics of half of the string length, /„ = Airn/l, as 



p = E^"' 



(1) 



n=l 



where P„ 



oc n 



Typically values of q are 4/3 for loops which generate cusps [Tg| , a point on the loop which travels at 
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FIG. 5: Pulsar constraints on Gfi/c as a function of dimensionless loop production size a. The various lines show limits from 

different observational constraints on t he g ravitational wave background 
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a — 60 G/i/c^ (long-dashed). 



2 < 2 X 10"" IgU, 
We also show the relation 



the speed of light at some point during its evolution, and 2 for a square loop which has 4 kinks on it [55|. However, this 
simplification has a significant drawback: it assumes that the string solution is not affected by radiation backreaction. 
This would mean that the cusps and kinks, both of which are points on the string moving at the speed of light, are 
not rounded off by the radiation they produce. 

The issue of gravitational radiation backreaction is difficult to address quantitatively, but it might be possible to 
glean some understanding of the process by studying Goldstone boson radiation from global strings instead. Although 
it is typically much stronger, this radiation is in many other ways similar to gravitational radiation, but in contrast it 
can be studied using a field theoretic approach. It has been shown that kinks on global strings are rapidly smoothed 
out by Goldstone boson radiation, and an initial kink perturbation is more or less indistinguishable from a sinusoidal 
perturbation after a few timesteps [63| ; a similar effect is expected for cusps. This suggests that the radiation spectrum 
from a loop will be modifed cither by a change in the power law, or strong cut-off at large n. 

In order to make contact with bounds on a stochastic background of gravitational waves from pulsar timing we 
need to estimate the spectral energy density of gravitational waves 



n„ 



f dPg 
Pg df 



emitted by the network. This can be related to the dimensionless strain by 

/100Hz\ 



1.3 X 10" 



n^K^ 



\ f J 
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In ref. [531 the effects of varying the power spectrum of radiation from the strings were studied. It was shown 
that the observed gravitational radiation background in the region relevant for pulsar observations (/ '^ 10"^ Hz) 
was sensitive spectrum of radiation emitted by the strings. In particular, it was shown that the amplitude of the 
background at these frequencies was around a factor of 3 higher when q ~ 4/3 compared to q = 2. A similar result 
was obtained by truncating the q = 4/3 spectrum at some point n* ~ 1000. The reason for this is that the emission 
at this particular frequency receives contributions from loops formed in the radiation era and also those from the 
matter era if they emit into high harmonics. Whatever one does with the spectrum of radiation emitted by loops from 
the matter era, the radiation era contribution will remain and, therefore, given the uncertainties in the spectrum of 
radiation from loops discussed above, we will take a conservative approach and just use the radiation era contribution 
as a lower bound on the signal for a given G^/c^. Tighter, but less reliable, constraints can be found by including 
the matter era contribution. 

We will assume that a is undetermined and use the analytic form for the radiation era spectrum presented in 
ref. [13 

^2 _ , ,, .. ,n-4GM fl-{vL,)/c'\ (1 + 1.4t)3/2 - 1 



n^h' = 1.17 X 10-*^ — ' '^^" ^ ' ■ ' ^ (4) 



where x = a(?/(TG^) and fim is the total matter density relative to the critical density. In what follows we use 
parameters measured from the Nambu simulation to give ^rad, (^rad) ^^^ ^^ ^*^t ^m ~ ^-^^ ^'^ ^^^^ ignore the 
uncertainties in these parameters since their cummulative effect is likely to be less than the uncertainties imposed by 
the lack of knowledge of a and the radiation spectrum. 

The constraints on G^ijc? are presented in Fig. [5] for F = 60 and as a function of a for a number of different 
possible values for the limit on fig/i^. Probably the most reliable of the published limit is ^^Y? < 2 x 10~* [64| . 
There is stronger published limit of 2 x 10~^ 1651 . although this is controversial [601. ^^ have also included the value 
of 9.3 X 10^* [66[ which is that used in ref. [531. ^^ each case the limit is flat for a ^ TG^/c^ and is oc a^^ for 
a ^ TGfi/c^; this asymptotic behaviour can be understood easily from Q. Using the limit of ref. [6^ we find that 
G/j/c^ < 7 X 10~'^ for ac^/{TGfi) < 1 and Gi-i/c^ < 5 x 10""/a for ac^/iVGii) > 1. In the former case this is 
weaker than the bound from the CMB and in the latter could be much stronger dependent on the value of a. 

V. CONCLUSIONS 

There are a wide range limits on Gfi/c^ presented in the literature, often using identical or at the very least similar 
data. In this paper we have tried to clarify these issues and present what we believe are current constaints from the 
CMB and large scale structure data, and also from pulsar timing. We have presented evidence that the string spectra 
which come from different simulation techniques can all be modelled using the USM with parameters measured in 
simulations and deduce limits of Gfi/c^ < 2.6 x 10~^ if the Nambu simulations are correct and Gn/c? < 6.4 x 10""^ if 
the AH simulations are correct. These use CMB+SDSS data and gstr as the parameter. We have also included data 
from BBN and used \og[G^/c^] as the parameter with similar but quantitatively different conclusions. 

The limits from pulsar timing could already be much stronger than those from the CMB, for example, if a is large. 
However, the uncertainties in interpretation of the pulsar limits on Q.gh? are much more severe. At this stage, it seems 
safest to use the low a limit of Vtgh? < 8 x 10"^, which is one that is unlikely to have to be taken back. 

As a final point we should note that a network of cosmic strings has another possible observational consequence 
in the CMB, albeit one of lower amplitude: the creation of B-mode polarization |67H63 l whose detection to very low 
values of G^ might be possible if lense cleaning is found to be possible [ZOl . This is due to the fact that a substantial 
fraction of the anisotropics created by strings are vector modes. In Fig. [S] we have presented the B-mode spectra 
for the Nambu and AH models with the value of G^l/c^ set to their repsective upper limits. These spectra are very 
different to those predicted in by infiationary models parameterized by the scalar-to-tensor ratio, r which are also 
included in the figure. In fact, the spectra for £ > 30 are similar to that produced by gravitational lensing particularly 
for tte Nambu case. This is because the B-mode polarization is created along the line of sight, albeit on different 
scales to the gravitational lensing pf the primary E-niode signal. These spectra comprise a white noise portion and 
a peak corresponding to the dominant scale which is on slightly larger scales in the AH model and, by coincidence, 
very similar to the lensing for the Nambu model. For the Nambu cas this similarity will make them difficult to detect 
unless one can detect the analogue of the reionization bump which peaks around (. ~ 10. It is, however, interesting 
since most infiationary models which produce strings create very low levels of gravitational waves. It has been shown 
that it is possible to discriminated them in the AH case [Tlj . 
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FIG. 6: _B— mode polarization power spectra from strings, tensors and gravitational lensing. For the string spectra, we use 
values of G^/c? corresponding to the 2a upper limit from CMB+SDSS data, that is. G^,/c? = 2.6 x 10~^ for the USM Nambu 
model (solid) and 6.4 x 10^^ for the AH case (dotted). We also show the inflationary primordial tensor spectrum with r = 0.1 
at fc = 0.05 Mpc~^ (short dashed) and r — 0.01 (long dashed). Finally, we show the gravitational lensing spectrum generated 
from _E— mode mixing (dot-dash) expected in the inflationary model. 
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